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interaction  gaps  occuring  at  two-level  and  three-level  resonances.  The  lowest  branch  lies  below 
the  bulk  TO  energy.  These  results  are  found  in  good  agreement  with  recent  experimental  data. 


I.  INTRODUCTION 


The  interface  and  confined  bulk  phonon  modes  in  heterostmctures  and  superlattices  have 
been  discussed  theoretically^  and  observed  experimentally^  in  the  past.  Very  recently,  optical 
phonon  modes  supported  by  a  semiconductor  double  heterostnicture  (DHS)  have  been  solved 
independently  in  the  long  wavelength  linut  by  two  groups.^’^  It  is  found  that  there  exist  two  types 
of  phonon  modes,  the  confined  bulk  modes  and  the  interface  modes.  The  former  modes  can 
be  either  longitudinal  optical  (LO)  or  transverse  optical  (TO)  with  frequencies  o^nd  (^iden¬ 
tical  to  those  of  the  bulk  excitations,  while  the  interface  (IN)  mcxles  may  have  frequencies  be¬ 
tween  (uf^and  (^  Eigenvectors  and  dispersion  relations  for  all  these  modes  can  be  found  in 
Refs.  3  and  4.  The  Hamiltonian  operator  describing  electron-optical-phonon  interactions  has 
also  been  derived  to  study  the  polaron  mobility  and  magnetophonon  resonance  spectra,^  and 
the  polaronic  states  in  a  DHS.^  The  importance  of  interface  modes  are  clearly  demonstrated 
in  these  calculations. 

Recent  experiments  of  far-infrared  photoconductivity  measurements^  carried  out  for  a 
donor  impurity  doped  at  the  center  of  GaAs  quantum  wells  in  a  GaAs/AlGaAs  multiple-quan¬ 
tum-well  (MQW)  structure  have  revealed  a  surprising  apparent  pinning  phenomenon.  The 
resonant  magnetopolaron  effect  on  the  ls-^+  transition  energy  is  found  to  be  around  on*,  in 


1 


contrast  to  the  bulk  case  which  has  been  well  understood  theoretically.  As  the  electron  does 
not  couple  to  TO  phonons,  the  data  have  been  regarded  as  a  significant  deviation  from 
reasonable  expectations  based  on  the  Frdhlich  model  of  interaction  with  the  bulk  GaAs  zone- 
center  LO  phonons,  both  m  magnitude  and  in  field  dependence.  The  possibility  of  zone  folding 
of  the  optical  phonon  dispersion  has  been  suggested^  to  explain  these  results.  Since  all  the  bulk 
modes  are  strictly  confined  in  heterostructures  or  superlattices,  there  cannot  be  any  traveling 
wave  in  the-z-direction  unless  the  adjacent  layers  in  the  MQW  structure  have  matching  optical 
phonon  frequencies. 

We  present,  in  this  article,  a  detailed  account  of  the  theory,  for  which  a  brief  outline  has 

o 

been  provided.  The  ls-2p  +  transition  energy  is  cai  julated  as  a  function  of  the  applied  magnetic 
field  for  a  magnetopolaron  bound  to  a  hydrogenic  impurity  in  the  quantum  well  of  a  DHS.  The 
electron-phonon  interaction  Hamiltonian  derived  in  Refs.  4  and  5  is  treated  as  a  perturbation 
on  tb:  hydrogen-like  impurity  atom  confined  in  a  quantum  well  under  strong  magnetic  fields. 
Since  the  unperturbed  system  cannot  be  solved  exactly,  the  variational  method  is  applied  by 
choosing  trial  wave  functions  with  proper  symmetry  properties.  It  is  found  that  interaction  gaps 
appear  in  the  transition  energy  which  then  breaks  into  separate  branches.  The  lowest  branch 
lies  below  the  TO  frequency  orr.  Vf^th  just  one  adjustable  parameter  in  the  trial  wave  functions, 
our  results  are  already  in  good  agreement  with  the  newly-improved  data. 


Very  recently,  the  experiment  has  been  improved  and  extended  in  a  series  of  measure¬ 
ments.'  Both  the  two-level  resonance  (the  ls-2p+  transition  energy  AE  is  equal  to  the  bulk  LO 
phonon  frequency  (ul)  and  the  three-level  resonance  (AE  is  equal  to^L  plus  the  energy  dif¬ 
ference  between  levels  2p.  and  Is)  are  observed  by  tuning  the  magnetic  field.  The  data  appear 
to  deviate  smoothly  from  the  calculated  transition  energy  in  the  absence  of  electron-phonon  in¬ 


teractions,  and  the  deviation  is  generally  smaller  than  what  is  reported  in  Ref.  In  addition, 
two  interaction  gaps  corresponding  to  the  two-level  and  three-level  resonances  are  observed. 


and  many  more  points  are  measured  above  the  first  gap.  As  is  shown  later,  these  data  agree 
quite  well  with  our  theoretical  curve,  where  possible  sources  of  small  corrections  are  given. 
We  outline  the  theory  in  Sec.  II  and  describe  the  method  of  calculation  in  Sec.  III.  Our 


results  are  discussed  along  with  conclusions  in  Sec.  IV. 


n.  THEORY 


Consider  a  donor  impurity  at  the  center  of  the  GaAs  quantum  well  of  width  d  in  a 
GaAsA^lGaAs  DHS  system.  A  magnetic  field  is  applied  along  the  growth  direction.  For  con¬ 
venience,  we  introduce  the  two-dimensional  vectors  ic  and  ^  such  that  k  «  (!c,q)  and  T  «  0o,z) 
for  the  phonon  momentum  and  electron  position,  respectively.  The  electron  momentum  is 
denoted  by  ke  »  (1^  |,  kz).  The  total  Hamiltonian  can  be  written  as 


where 


(2a) 


is  the  Hamiltonian  for  a  hydrogenic  impurity  confined  m  the  square  well  Vb(z)  given  by 


|z|  <  d/2 

\z\  >  d/2. 


(2b) 


The  free  phonon  Hamiltonian  is 
Hph  =  Hlo  +  Hin 


Hlo  =  ^ 

Z,m  ^ 


(3a) 

(3b) 

(3c) 


In  Eqs.  (3),  we  have  introduced  the  index  v  to  label  the  material,  with  v  =>  1  for  the  well  and  v 
«  2  for  the  barrier.  We  have  also  defined  the  creation  (annihilation)  operators  afn^(K)  (am(A:)) 
for  the  confined  modes  and  a^(j,«j(/c)  (aij.aj(ic))  for  the  symmetric  and  antisymmetric  interface 
modes  of  phonons.  They  obey  the  commutation  relations 

[ao(;0.aj(^)]  =  «a4«(«-/?) 


(al(«),a{(«')]  =  (a«(ic),a^(;?)]  =  0. 


(4a) 

(4b) 


The  mteraction  Hamiltonian  He-ph  is  taken  directly  from  Ref.  S.  It  consists  of  two  terms: 
the  electron  interaction  with  confined  LO  modes,  and  the  interaction  with  interface  modes.  As 
has  been  shown  in  Ref.''  3,  5,  contributions  of  lattice  vibrations  outside  the  well  to  the  polaron 
effect  are  significant  only  when  the  well  width  d  is  extremely  small.  For  d  2  100  A  they  are  com* 
pletely  negli^ble.  Therefore  we  have 


\diere 


(5a) 


(5b) 


represents  the  electron  interaction  with  confined  bulk  LO  phonon  modes  m  the  quantum  well 
in  which  q  nui/d  is  quantized,  and 


(5c) 


represents  the  electron  interaction  energy  with  interface  phonon  modes  in  the  well.  The  nor¬ 
malization  constants  are  given  by 


Ad  +  (mv/dy  \eooi  ec^^  J 


?rc- 


Ak  i\  tanhjKd/2) 

IB 

'  i4«  ?i  COth(Kd/2)  +  £2  ’ 


^  . -i  i 

•(6b) 
(6c) 


vt^ere  A  stands  for  the  interface  area,  eoov  and  Cov  denote  the  optic  and  dielectric  constants  of 
material  v,  and  ev(£o)  is  defmed  by 


1  _  1 _ 1  I 

iy(u/) .  ri,(w)  —  £j„  ^t'(^)  “ 


with  the  dielectric  function  of  material  v  given  by 
€^{u)  =  -  w^)/(wt„  -  w’).  i 


OD  OF  CALCULATION 


(7) 

(8) 


We  now  proceed  to  calculate  the  transition  energy  by  perturbation  theory.  For  the  unper¬ 
turbed  system,  we  take  the  hydrogenic  impurity  m  the  well  with  a  magnetic  field  applied  normal 
to  the  interfaces.  Thus  Ho  >■  He  +  Hph  and  He-pii  is  treated  as  a  perturbation.  The  energy 
levels  are  calculated  by  means  of  Wigner-Brillouin  perturbation  theory  to  second  order.  The 
result  is 


e<(F)  =  f:.(B)-l- 


^  Y^j2,r  I  <  >  P 

i2x)^^  uiB)-E„{B)^hu,{Ky 


(9) 


whCTe  ]fe(B)  stands  for  the  corresponding  unpmurbed  energy,  and  the  state  |n>  |n/;  N> 

is  specified  by  the  atomic  state  ol  and  the  phonon  number  N. 

It  is  noted,  however,  that  even  the  unperturbed  problem  is  not  exactly  solvable.  Thus  we 


employ  the  variational  method  to  determine  the  unperturbed  energy  levels.  The  trial  wave  func¬ 
tion  is  taken  to  be 


rp  =  f{z)G{p,z,4>).  1 

The  function  f(z)  in  Eq.  (9)  is  the  solution  to  the  square-well  problem.  It  is  given  by 


(10) 


f(A  -  !  \z\  >  d/2 

^  \  >lexp(-I:;iz|),  izi>d/2 


where  the  wave  numbers  kz  and  kz’  are  related  to  the  first  electron  subband  energy  Ei  by 

k,  =  yj2m,Eilh'^,  k'^  =  ^J2m,[y^  -  '  (12) 

The  function  is  chosen  to  reflect  the  symmetry  properties  of  the  system.  It  is  given  by 

G{p,  z,  0)  =  p’"*'  exp(im.^  -  yp^/A  -  yi^z^/4),  i  (13) 


I 


1 


where  we  have  defined  the  parameter 

y-huc/2Ry,  I  (  (14) 

with  the  cyclotron  frequency  cuc  “  eB/meC  and  the  effective  Rydberg  Ry  =  mcC*/2k^  for  the  im- 
■pirity!~The  quantum  numbeTln  specifies  the  impurity  levels  such  that  »  1,0, -1  for  the  2p+, 
Is,  2p.  levels,  respectively.  The  variational  parameter  { is  determmed  by  minimizing  the  energy 
of  the  level  in  question.  All  the  other  levels,  for  the  particular  samples  used  in  these  measure¬ 
ments,^  can  be  neglected  in  these  experiments.  Hence  the  Hilbert  space  in  our  perturbation  cal¬ 
culation  is  truncated  to  only  the  three  states  Is,  2p.  and  2p+. 

As  the  matrix  element  is  generally  small  compared  to  the  transition  energy,  the  perturba¬ 
tion  energy  in  (9)  becomes  appreciable  only  when  the  energy  difference  in  the  denominator 
matches  the  phonon  energy.  This  implies  immediately  that  the  electron-phonon  coupling  has 
negligible  influence  on  the  Is  energy  level.  Hence  it  is  suffldent  to  calculate  AE  =  ezp-t-  -  Els 
for  the  transition  energy.  Furthermore,  we  note  thitf  significant  contribution  from  the  second 
term  in^to)  is  expected  around  AE  w  Ezp.  •  Eu  +  as  well  as  around  AE  » 

As  can  be  seen  from  Eqs.  (^,  the  calculation  of  energy  levels  involves  the  evaluatiuon  of 
four  matrix  elements 


M\(l)  -<  2p+,0|/fjt_i,*lla,l  > 


A^sCt)  =<  2p+,01/f*«/Af|l#,  1  > 
Mi{y)  =<  2p+,01/f,_i,o|2p-,  1  > 

^4(7)  =<  2p+,0|ff,_//vl2p>,l  >  . 


y 


1.:^  di. 
'I 


(15) 


Explicit  expressions  for  these  matrix  elements  are  worked  out  in  the  Appendix.  Here  we  just 
present  the  results  as  a  function  of  the  applied  field  in  Fig.  1.  It  is  observed  that  in  both  the 
two-level  and  three-level  resonance  cases,  the  interface  modes  always  yield  larger  matrix  ele¬ 
ments  than  the  corresponding  confined  bulk  modes,  and  that  Mt  is  the  largest  and  is  generally 
3  to  4  times  larger  than  the  others.  This  is  understandable  because  the  first  Bohr  radius  of  the 
impurity  is  of  the  order  of  the  well  width  and  the  p-state  wave  functions  extend  out  much  farther 
m  the  xy-plane  than  the  s-state  wave  function. 


FIG.  1.  Absolute  value  of  the  matrix  elements  of  electron-phonon  interactions  calculated  as  functions  of  the 
applied  magnetic  field  in  the  GaAs  well  of  a  GaAs/GaAlAs  double  heterostructure. 


The  energy  of  2p+  level  is  then  given  by 

I 

<2p+(7)  =  ^2p+(7)  +  T{-r),  j  (16) 

where  T(y)  *  Ti  +  T2  +  T3  -h  T4  represents  the  perturbation  energy.  More  explicitly,  we  have 


7i  = 


2a{hui,) 
rK, 


it-/ 


I  <  2p^.|e'^  ^sin[mT(*  -I-  i)]jla  >  |» 


U2Sl(7) 
f2p^  —  Eu  —  fML 


r 


dt 


f* -I- (miro,/d)* 


(16a) 


2 


^  r  hu,j(K)  gj(K)  I  <  2p+|e*^  ^cosh(/cz)|l5  > 


cosh^(Y)  «2p+  -  El,  -  hu,j{K) 


hu)aj{K)  gj(/c)  I  <  2p.t.|e*^^sinh(/cz)|l 
K  sinh^(^)  C2p+  -  Eu- 

?v  hu..At)fiO 


cosh^{id/aa) 


2p.t.|e*^^sinh(/cz)|l5  >  P  '1 
E2p+  -  Eu-  hLJaj{K)  } 

,^U2p+ -  ^^1,  - 


f  ,2  I  <  2p+|e‘^-^sin[mff(^  -  f)]|2p-  >  ^ 

aihutY  ^  ie^MSlh)  r  t(2y-i^e-*^f^y)^ 

Ivk-d  2^  _  fc,.,,  /  ¥2 


r  /;  4 


25*pd  ^  f2p+  ^  E2p_  -  huL  Vo  -f  (mvae/tt^ 

—  [  f  I  2p.f[e*^'^cosh(<cz)|2p-  >  p 

''^^0  J  C0Sh^(£//2)  ^2p+  —  E2p^  —  hlji,j{K) 

^aj{K)  gj{K)  I  <  2p+|e'^  ^sinh(Kz)|2p-  > 

«  sinh^(d/2)  f2p+ -  ^^2p_  -  ^«2(k)  / 


=  ^  Y^r  hu.i(t)Mt,y)  ] 

8*'* /o  7cosh^(<d/ao)  ^[<2p>  —  ^2p_  —  !  ’ 

where  and  I2  are  the  variational  parameters  for  Is  and  2p.  states,  respectively,  and  are  deter¬ 
mined  by  minimizing  the  energy  level  in  question,  ao  is  the  first  Bohr  radius  of  the  impurity  atom, 
and  Kp*^  =  jJi/2me<UL)^^  measures  the  size  of  the  polaron  due  to  LX)  phonons.  We  have  also 
defined  in  Eqs.  (16)  the  variable  t  =  Kao  and  the  functions 


5m(7)=  /  dze-T'(«?+«>^"*/^sin[rJ^z  + ~  co8^(fc,z)  (17) 

J-d/2a,  2 

7)  =  9j(t)  <i2:e-'>'(«? +€?)**/<  co8h(tz)  cos^(ir,z)j 

(0  =  f*i/[€i(wy)  tanh(td/2a,)  +  «2(w/)]- 

The  integrals  in  (16)  are  still  very  difficult  to  evaluate.  We  adopt  a  numerical  iteration  pro¬ 
cedure  to  compute  the  transition  energy.  Let  X(y)  =  £2pA  *  and  r(y)  »  E2p^-  Eu.  Then 
Eq.  (^  can  be  rewritten  as  /  '  •  - 

X{y)  =  r(7)  +  r(7,  X).  (20) 

The  function  r(y)  represents  the  transition  energy  in  the  absence  of  the  electron-phonon  inter- 


actions  and  is  already  known  from  the  variational  calculations.  Therefore,  we  start  the  iteration 
with  X  =  r(y)  in  T  for  a  fixed  y.  A  new  X  value  is  obtained  from  (20)  and  employed  to  compute 
a  new  T.  The  procedure  goes  on  until  self-consistency  is  achieved.  The  iteration  procedure 
repeats  when  y  changes,  and  eventually  the  transition  energy  is  found  as  a  function  of  y,  which 
is  proportional  to  B. 


IV.  RESULTS  AND  DISCUSSION 


Results  of  our  numerical  computation  are  presented  in  Fig.  2  in  which  the  ls-2p+  transi¬ 
tion  energy  is  plotted  as  a  function  of  B.  The  parameters  used  in  the  numerical  work  are  d  = 
125  A,  me  =  0.067  mo,  a©  =  98.7  A,  Ry  =  5.83  meV  and  V©  =  2,30  meV.  The  solid  line  indudes 


FIG.  2.  Comparison  of  experimental  data  with  the  calculated  Is-Zp^'  transition  energf  as  a  function  of  the 
applied  field  B.  The  solid  line  represents  results  from  the  present  theory  including  nonparablk  band 
mass.  The  dashed-dotted  line  does  not  include  the  nonparabolidty  correction,  and  the  dashed  line 
b  the  transition  energy  in  the  absence  of  electron-phonon  interactions  calculated  from  Gieene-Bajaj 
wave  functions  with  nonparabolic  mass  corrections  induded. 


reseats  the  transition  energy  without  electron-phonon  interactions.  It  is  calculated  from 
Greene-Bajaj  wave  functions^^  and  corrected  for  the  nonparabolicity  effect.  Experimental  data 
are  taken  from  a  series  of  measurements  reported  in  Ref.  9.  It  is  striking  that  the  theoretical 
curve  breaks  into  three  branches  separated  by  two  gaps,  in  good  agreement  with  experimental 
data.  The  pinning  effect,  or  the  unperturbed  energy  level  repelling,  is  apparently  a  resulti^f 
the  strong  resonance  interaction  when  the  denominator  of  the  second-order  perturbation  energy 
vanishes.  As  the  electron  does  not  interact  with  TO  phonons,  the  appearance  of  Jl^^aps  can 
only  be  attributed  to  the  coupling  of  the  electron  with  interface  phonons  which  oscillate  at  fre¬ 
quencies  between  the  bulk  cut  and  cul. 

Further  studies,  both  theoretical  and  experimental,  are  necessary  for  a  complete  under¬ 
standing  of  the  problem.  The  determination  of  the  transition  frequency  from  experimental  data 
is  difficult,  as  the  observed  intensity  distribution  in  the  resonance  region  deviates  greatly  from 
the  Lorentz  shape.  Theoretical  study  of  the  influence  of  electron-phonon  interactions  on  the 
transition  probability  is  needed  for  detailed  comparison  with  the  experimental  line  shapes  and 
observed  rapid  decrease  in  intensity  in  this  energy  region.  Furthermore,  as  we  have  noted  pre¬ 
viously,  the  trial  function  (13)  has  only  one  adjustable  parameter.  More  accurate  functions  must 
be  constructed,  for  more  accurate  results. 

In  conclusion,  we  assert  that  it  is  the  interaction  of  the  electron  with  interface  phonon 
modes  that  change  the  polaronic  properties  in  the  reduced  geometry.  These  modes  should  al¬ 
ways  be  included  in  the  interpretation  of  such  experiments.  Since  the  coupling  constant  of  the 
bulk  Frbhlich  type  cannot  be  clearly  defmed  for  the  interface  modes,  polaronic  phenomena  ob¬ 
served  in  the  reduced  dimensionality  should  Mbe  analyzed  by  varying  the  coupling  strength 


in  terms  of  the  bulk  Frdhlich  interaction. 
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The  integration  of  MD  matrix  elements  can  be  evaluated  in  the  following  manner.  We  start 
with  Ml  »  Iffl  Mira  where 


M\m  =<  2p+le’*  ^sin 


,z  1 . 


|ls>. 


(Al) 


Using  the  wave  functions  from  Eqs.70)-(l^w^avc 


Mira  =  clyJ^SM  -  t  sin  ,  (A2) 


where  we  have  defined  the  function 

rm 


(7)=  / 

J-rf/2 


t 


-"*(—+ d 


(A3) 


^  ' 


As  we  have  already  noted  before,  the  integration  of  z  can  only  be  done  numerically.  Here,  we 
shall  work  on 


Jo 

«/nsi 


f  d(p{cos 
Jo 


0  2(n  +  l)!!  JO 

~  (-l)"(a,/cp)^"^->  [2(n-H) -!]«!, 


(2n+l)!  2»+i(n  +  l)  ^  ^ 

-  ■’  L  1 1)!-!  = 


(A4) 


where  Ji(aoKp)  is  the  Bessel  function.  The  sine  term  in  (A2)  vanishes  by  direct  integration. 
To  evaluate  |  Mim  I  ^  we  need 


If 


^2 


(A5) 


<ipp’e-^'*/Vi(a,/cp) 

0 

(n+l)!n!  J, 

=  v'-tdii-M  ”*‘i  fii 
2  J  2I2J 


A  combination  of  (A2),  (A3)  and  (A5)  then  yields  the  last  line  of  Eq.  (16a). 

In  the  calculation  of  M2,  we  need  only  to  consider  the  symmetric  term,  as  the  antisymmetric 
term  vanishes.  Thus 


A/j  =<  2p+le'*'^co8h(K*)|l«  >, 


(A6) 


vdiich  has  the  same  form  as  Mim  except  for  the  z-dependence.  Hence,  the  above  procedure  ap¬ 
plies  provided  that  SmO*)  w  replaced  by  the  function  fj(t,y)  as  given  by  (18). 

The  matrix  element  M3  «  £  Mam  where 
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